#include <iostream>

#include "Integral.h"
#include "Density.h"
#include "Gradient.h"
#include "voronoi.h"

int main() {
    peopleDensity = new double[N * N];
    industryDensity = new double[N * N];
    getArray(getPeopleDensity, &peopleDensity);
    getArray(getIndustryDensity, &industryDensity);

    InitData();

    DataList p;
    for (int i = 0; i < K; ++i) {
        p.push_back(0.1 + i / (2.0 * K));
        p.push_back(0.2 + i / (2.0 * K));
    }
    
    double sigma = 0.000001;
    double epsilon = 0.000001;

    GradSearch(p, sigma, epsilon);

    delete [] peopleDensity;
    delete [] industryDensity;
    delete [] vertices;
    return 0;
}
